Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 11/12/01 



THE MOMENT OF INERTIA OF THE BINARY PULSAR J0737-3039A: CONSTRAINING THE 

NUCLEAR EQUATION OF STATE 

I. A. Morrison 1 , T. W. Baumgarte 1 ' 2,3 , S. L. Shapiro 2 ' 4 and V. R. Pandharipande 2 

Draft version February 2, 2008 

ABSTRACT 

We construct numerical models of the newly discovered binary pulsar J0737-3039A, both with a fully 
relativistic, uniformly rotating, equilibrium code that handles arbitrary spins and in the relativistic, 
slow-rotation approximation. We compare results for a representative sample of viable nuclear equations 
of state (EOS) that span three, qualitatively different, classes of models for the description of nuclear 
matter. A future dynamical measurement of the neutron star's moment of inertia from pulsar timing 
data will impose significant constraints on the nuclear EOS. Even a moderately accurate measurement 
(< 10%) may be able to rule out some of these competing classes. Using the measured mass, spin and 
moment of inertia to identify the optimal model computed from different EOSs, one can determine the 
pulsar's radius. 

Subject headings: Gravitation — stars: neutron — pulsars: individual (J0737-3039) — stars: rotation 



1. INTRODUCTION 

The recently discovered binary neutron star system 
J0737-3039 has many very interesting and exciting proper- 
ties (e.g. Burgay et. al 2003; Lyne et. al 2004). Both stars 
in J0737-3039 have been identified as pulsars, the binary 
is about ten times closer to the Earth than the Hulse- 
Taylor binary B1913+16, and the binary separation is sig- 
nificantly less than that of B1913+16 (Burgay et. al 2003). 
As a consequence, relativistic effects affecting the orbital 
dynamics play a much larger role in J0737-3039 than in 
all previously discovered binary neutron stars. One such 
effect that may become measurable soon from pulsar tim- 
ing data is spin-orbit coupling (Damour & Schafer 1988; 
Wex 1995). In particular, Lens-Thirring precession will 
yield the moment of inertia of the more rapidly spinning 
pulsar J0737-3039A (Lyne et. al 2004; O'Connell 2004). 
Such a determination would constitute the first time that 
the moment of inertia of a neutron star (or any other 
star) has been measured by dynamical means. The re- 
cently launched Gravity Probe B experiment will measure 
Lens-Thirring as well as geodetic precession of a gyroscope 
in orbit about the Earth. Geodetic precession has been 
measured in the binary pulsar system B1534+12 (Stairs, 
Thorsett & Arzoumanian 2004). 

Since the star's gravitational mass and spin frequency 
are well determined, a measurement of the moment of in- 
ertia could provide a stringent constraint on the nuclear 
equation of state (EOS). Using the measured mass, spin 
and moment of inertia to select the optimal candidate from 
a suite of theoretical models constructed from competing 
EOSs, one can also determine the pulsar's radius. 

In this paper we point out that even a measurement of 
the moment of inertia of modest accuracy (< 10%) may 
help differentiate between competing approaches that are 
presently adopted for the description of nuclear matter and 
the construction of nuclear EOSs (for a general review of 



modern EOSs see Lattimer & Prakash 2000, 2001). In Sec- 
tion 2 we briefly review three qualitatively distinct classes 
of EOSs; in Section 3 we adopt representative EOSs in each 
of these classes to construct numerical models of J0737- 
3039A and determine its moment of inertia; in Section 4 
we briefly discuss our findings. 

2. EQUATIONS OF STATE 

In this paper we consider six different nuclear EOSs. 
They are chosen to represent three different classes of 
EOSs, distinguished by the formalism used in their con- 
struction (see discussion below). All six of these EOSs give 
maximum gravitational masses that exceed the mass of the 
most massive neutron star (PSR J1518+4904, M/M = 
1.56; Thorsett & Chakrabarty 1999) residing in the six 
known binary radio pulsar systems. These systems provide 
neutron stars with the most accurately determined masses. 
Some pulsar masses determined with greater uncertainty 
(e.g., PSR J0751+1807, believed to have a low-mass white 
dwarf companion and an estimated mass M/Mq ~ 2.2; 
see, e.g., Table 1 in Stairs 2004) are greater than the max- 
imum allowed masses for some of the chosen six EOSs (see 
Table 1 below). If confirmed, these objects would rule out 
these particular EOSs, or require adjustment of some of 
the parameters on which they are based. 

The first class of EOSs that we consider are results 
of nonrelativistic many-body calculations using realis- 
tic models of two- and three-nucleon interactions. The 
two-nuclon interactions are obtained by fitting observed 
nucleon-nucleon scattering data, and the three-nucleon in- 
teractions are based on the energies of light nuclei and/or 
properties of symmeric nuclear matter (SNM) with com- 
parable number of neutrons and protons. To represent 
these class I EOSs we choose two examples. The first is 
APR (Akmal, Pandharipande, & Ravenhall 1998), which 
adopts a modern two-nuclon interaction (A18), together 
with its relativistic boost correction and the UIX three- 
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nucleon potential. The second is FPS (Lorenz, Ravenhall, 
& Pethick 1993), which is based on an older two- nucleon 
interaction (U14), and phenomenological density depen- 
dent terms to represent the effects of three-nucleon forces. 
The main difference between these two EOSs is that FPS 
exactly reproduces the empirical properties of SNM, but 
not light nuclei, while the A18+UIX forces used by APR 
reproduce both the binding energies of light nuclei and 
SNM properties fairly accurately. 

An important property of the EOSs belonging to class 
I is that the pressure of pure neutron matter (PNM) at 
desities in the 10 11 — 10 14 g cm" 3 range is close to half 
of that of noninteracting neutron gas (Carlson, Morales, 
Pandharipande, & Ravenhall 2003). This follows from the 
known, large neutron-neutron scattering length, a <~ 18 
fm, which makes the dimensionless parameter akp, where 
kp is the Fermi momentum of PNM, larger than unity, 
and results in relatively smaller crusts, radii and moments 
of inertia for neutron stars with M < l.6M@. 

The second, class II, of EOSs are based on relativis- 
tic mean- field approximations. A mean-field Lagrangian 
is chosen to reproduce the empirical properties of SNM, 
and it is used to calculate the EOS of neutron star matter, 
which is closer to PNM, using the mean-field approxima- 
tion. SNM is charged, and in absence of Coulomb forces it 
is believed to be bound with an energy of ~ —16 MeV 
per nucleon, and to have zero pressure at a density of 
~ 2.5 x 10 14 g cm~ 3 . At lower densities its pressure is 
negative and it becomes unstable at densities less than 
~ 1.5 x 10 14 g cm~ 3 . On the other hand PNM has possi- 
tive energy and pressure at all densities as far as we know. 
Thus it is most challenging to extrapolate the properties 
of SNM to those of PNM. The relativistic mean field La- 
grangian contains no information on neutron-neutron scat- 
tering, such as the large value of the scattering length a, 
and therefore the pressure of class II EOSs approaches 
that of noninteracting neutron gas at desities less than 
10 14 g cm~ 3 . They yield neutron stars with relatively 
larger crusts, radii and moments of inertia. We choose 
MSI (Mullcr and Serot 1996) and GM3 (Glcndcnning & 
Moszkowski 1991) to represent class II. 

The class III of EOSs we consider assume that neutron 
stars are entirely made up of strange quark matter. This 
hypothetical form of matter is electrically neutral and has 
zero pressure at densities of order 10 14 g cm -3 . There is 
no information or observational evidence for this matter. 
We consider two examples: SQM1 and SQM3 (Prakash, 
Cooke, & Lattimer 1995). 

Tables of the comoving total mass density p, the rest 
mass density po & n d the pressure P of the chosen EOSs 
were kindly provided by M. Prakash and D. G. Raven- 
hall. To improve the convergence of our numerical codes 
we found it helpful to increase the number of data points 
using third-order interpolation between the raw data. 

3. THE MOMENT OF INERTIA 

Adopting the EOSs discussed in Section 2 we compute 
the moment of inertia of the pulsar J0737-3039A by per- 
forming two independent relativistic calculations. We as- 
sume that the pulsar is rotating uniformly; while nascent 
neutron stars can form with appreciable differential ro- 
tation, viscosity and magnetic fields will combine to drive 



them to uniform rotation on a timescale much shorter than 
the age of the pulsar inferred from its spindown, 210 Myr. 
In addition, the precision of the underlying pulsar clock 
mechanism strongly argues in favor of solid-body rotation. 

In the first calculation, we adopt the numerical code 
of Cook, Shapiro & Teukolsky 1994, hereafter CST, to 
construct stationary models of uniformly rotating neutron 
stars. This code treats arbitrary spin rates without ap- 
proximation. We use it to compute models that have the 
gravitational mass M — 1.337Mq and the angular velocity 
£1 = 276. 8s~ 4 (corresponding to a spin period P = 27.7 
ms) of J0737-3039A. We refer to CST for details of the 
code, including how the code utilizes tabulated EOSs. The 
code computes the angular momentum J 



J 



(1) 



where T^ v is the stress-energy tensor for a perfect fluid, 

= ( P + P)uV +Pg» v , (2) 
and £^ is the Killing vector in the azimuthal direction 
reflecting axisymmetry. In (2) is the fluid four-velocity, 
which for uniformly rotating, axisymmetric stars takes the 
form u M = u*(l, 0, 0, O) where f2 = vfi ju 1 is the stellar an- 
gular velocity. The definition (1) then reduces to 



J = J(p + p)u t (g <l>4> u' l > + g c f >t u t )y/^jd i x 



(3) 



where g = det(g Q| g). From the angular momentum J the 
moment of inertia is determined as / = J/tt. 

A simple estimate shows that the stellar frequency fl = 
276.8 Hz is much less than the Kepler frequency at the 
stellar equator, fi max w (GM/i? 3 ) 1 / 2 . For all the EOSs 
we consider here f2/f2 max < 0.035. This implies that devi- 
ations from sphericity are very small, so that the moment 
of inertia can also be approximated from spherical stel- 
lar models. Our second calculation therefore adopts this 
slow-rotation approximation (Hartle 1967; Chandrasekhar 
& Miller 1974; Glendenning & Weber 1992; Ravenhall & 
Pethick 1994) which we briefly review as follows. 

In the slow-rotation limit in spherical polar coordinates 
the metric can be written as (G = c = 1) 



ds* = -e 2 *dt z + ( 1 - 



dr" 



- 2uor 2 sin 2 6dtd<j) + r 2 (d8 2 + sin 2 6d(f> 2 ) (4) 
where m(r) is the total mass within radius r 

— — = Anr z p (5) 
dr 

and where u>(r) = d<f)/dt\ ZAMO is the Lense-Thirring angu- 
lar velocity of a zero-angular momentum observer (ZAMO) 
with = 0. Up to first order in uj all metric potentials are 
spherically symmetric and only depend on the areal radius 
r. In the stellar interior, Einstein's equations reduce to 
d$ _ m + (4irr 3 p) 
dr 



1 - 2m/r 



(r < R), 



(6) 



and 



1 d ( , . duj \ dj r . . 



where following convention we define u> = — u> and 

2^ 1/2 
r 



1 



(7) 



(8) 
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Table 1 

Numerical results for J0737-3039A (M = 1.337M Q and Q = 276.8s -1 ) in the slow-rotation approximation 



Class 


EOS 


M max /M a 


J b 


I c 


I /MR 2 


o d 


M /M Q c 


ms/m q 1 


R s 


GM/c 2 R 


I 


APR 11 


2.20 


3.42 


1.24 


0.348 


9.61 


1.47 


1.37 


11.56 


0.1707 




FPS 1 


1.80 


3.14 


1.14 


0.358 


12.1 


1.48 


1.37 


10.92 


0.1809 


II 


MSI J 


1.82 


4.61 


1.66 


0.326 


6.81 


1.45 


1.34 


13.87 


0.1424 




GM3 k 


1.56 


4.14 


1.49 


0.325 


9.10 


1.46 


1.35 


13.12 


0.1502 


III 


SQMl 1 


1.56 


2.72 


0.982 


0.474 


15.1 


1.55 


1.44 


8.830 


0.2236 




SQM3 1 


1.99 


3.91 


1.41 


0.459 


7.00 


1.50 


1.40 


10.76 


0.1836 



a maximum allowed gravitational mass for non-rotating star 
b total angular momentum (10 47 g cm 2 s — eqn. (3) 
c moment of inertia (10 45 g cm 2 ), eqn. (14) 
d central mass-energy density (10 14 g cm -3 ) 
e rest mass 

f rest mass of companion, J0737-3039B (M = 1.25M , Cl = 2.27s- 1 ) 
3 circumferential radius (km) 

h A18 +Sv + UIX*; Akmal, Pandharipande, & Ravenhall (1998) 

1 UV14 + TNI; Lorenz, Ravenhall, & Pethick (1993) 

3 Field Theoretical; Glendcnnning & Moszkowski (1991) 

k Field Theoretical; Mullcr & Scrot (1996) 

1 Quark Matter; Prakash, Cooke, & Lattimer (1995) 



The stellar structure is determined by the Tolman- 
Oppenheimer-Volkoff eq. (Arnett & Bowers 1976) 5 

dP , ,m + 47rr 3 P 

— = -{p + P) . (9) 

dr 1 — zm/r 

Exterior to the stellar surface at r — R, defined as the 

radius R at which the pressure first vanishes P(R) = 0, 

the metric reduces to 

e 2 * =(l-—) (r>R), (10) 



r I 



and 



2J 



LU = 



(r > R) 



(11) 



where M — m(R) is the total gravitational mass. 

On the surface of the star the interior and exterior solu- 
tions are matched by satisfying the boundary conditions 

*{R) = n-~, (12) 



3 dr 



and 



1 



= - ln(l - 2M/R). 



(13) 



In practice this is accomplished by choosing arbitrary ini- 
tial values <& c and lo c at r — for the integration of (6) 
and (7) to the surface r = R. The boundary conditions 
are then met by multiplying the interior solutions by ap- 
propriate constants. 

The moment of inertia I = J/Cl can then be computed 
from eq. (3). Using ft — /u l and keeping only terms up 
to first order in ft and ui we find 



7^ 



(p + p)e * w 



■ 4 dr. 



(14) 



In this framework, the moment of inertia can be com- 
puted by solving ordinary differential equations, namely 
the Tolman-Oppcnhcimcr-Volkoff eqs. (5) and (9) together 
with eqs. (6) and (7), followed by a numerical quadrature, 
(14). This integration can be carried out with almost arbi- 
trary accuracy, so that that the error is dominated by the 
slow-rotation approximation. From (f2/^max) 2 ~ 10~ 3 we 
expect this error to be in the order of about 0.1 % (which 
we have verified by comparison with our calculations us- 
ing the CST code). For most EOSs our datasets are fairly 
sparse, which leads to an interpolation error that is closer 
to about 1 % for some EOSs, far larger than the intrinsic 
error arising from the slow-rotation approximation. We 
summarize our numerical results in Table 1, constructed 
using this approximation. All results obtained with the 
CST code were within 1 % of the results in Table 1, the dif- 
ference dominated by interpolation error within the EOS 
datasets, which is carried out differently in the two calcu- 
lations. (We found closer agreement for those EOSs for 
which more datapoints are available). We also point out 
that the CST code calculations showed that the oblateness 
of the stars is very small for all EOSs R e /R p — 1 < 10~ 3 , 
where R e and R p are the equatorial and polar radius. 



3 7 (1 - 2m/r)V2 ft 
This slow-rotation approximation for the moment of iner- 
tia neglects deviations from sphericity and is independent 
of the angular velocity ft. 

5 Note that in eq. (3.11) of Arnett & Bowers this expression has an incorrect exponent 



4. DISCUSSION 

It is evident from Table 1 that the moment of inertia I 
for a neutron star like J0737-3039A depends quite strongly 
on the EOS. Hence a measurement of I could impose new 
constraints on the EOS, refining some of the parameters 
that enter the models of nuclear matter, or even ruling out 
whole classes of models that currently exist. 
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Within each Class I and II, the predictions for the neu- 
tron star radius and moment of inertia are fairly simi- 
lar, but between these two classes, the values are quite 
different. Class I predicts a radius < 12 km and a mo- 
ment of inertia < 1.2 x 10 45 g cm 2 , while Group II pre- 
dicts a radius of about > 13 km and a moment of inertia 
> 1.5 x 10 45 g cm 2 . This implies that an observational 
measurement of / even with a modest accuracy (say 10 % 
or so) may distinguish these different groups, establishing 
whether "realistic" many-body nuclear models based on 
nuclear matter and experimental data (Class I) or mean- 
field models (Class II) provide better descriptions of neu- 
tron star matter. 

There is a much larger spread in our results for Class 
III, which seem to depend quite strongly on various free 
parameters in the EOS. These parameters affect the (sur- 
face) density at which the pressure vanishes, and hence 
the radius and moment of inertia of the stellar models. A 
measurement of / may serve to constrain these parameters. 

As an aside we point out that the pulsar's companion, 
J0737-3039B, which is also a pulsar, has a surprisingly 
small gravitational mass of only M = 1.25M Q . The cor- 
responding baryonic rest mass depends on the EOS and 
ranges from 1.34 M Q to 1.37 M Q for Classes I and II (see 



Table 1; also Fig. 3 in Morrison, Baumgarte & Shapiro 
2004). Most supernova simulations lead to remnants that 
are more massive than this (e.g. Woosley, Heger & Weaver 
2002), especially if general relativistic effects are included 
(e.g. Fryer 1999). It is likely that J0737-3039B originated 
from a close helium star-neutron star binary (Devi & van 
den Heuvel 2004; Willems & Kalogera 2004), which will 
clearly affect the dynamics of the supernova explosion and 
hence the remnant mass. It is therefore too early to draw 
conclusions, but any future simulations of supernovae in 
close binaries ought to be able to produce such small mass 
remnants. That outcome may be possible for some EOSs 
but not for others, which would put a further constraint 
on the nuclear EOS. 
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